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Abstract 

Representations for the generating functionals of static correlators of z-components 
of spins in the XY and XX Heisenberg spin chains are obtained in the form of sums 
of the fermionic functional integrals. The peculiarity of the functional integrals in 
question is because of the fact that the integration variables depend on the imaginary 
time "automorphically" . In other words, the integration variables are multiplied 
with a certain complex number when the imaginary time is shifted by a period. 
Therefore, the corresponding boundary conditions at the ends of the imaginary 
time segment are not of the form corresponding to fermionic, or bosonic, variables 
taken in the Matsubara representation at nonzero temperature. In fact, one part of 
sites of the models corresponds to the integration variables which are subjected to 
the unusual boundary conditions, while the variables on the other sites depend on 
the imaginary time conventionally, i.e., as fermions (or bosons). Thus a situation, 
when an "automorphic" boundary condition is the same for all sites of a chain spin 
model, is generalized. The results of the functional integration are obtained in the 
form of determinants of the matrix operators which are regularized by means of 
the generalized zeta-function approach. The partition functions of the models and 
certain correlation functions at nonzero temperature are obtained explicitly thus 
demonstrating correctness of the functional integral representations proposed. 
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1 Introduction 



The correlation functions of quantum models that are solvable via the Bethe ansatz 
method [TJ can be represented, in the thermodynamic limit, as the Fredholm determi- 
nants of certain linear integral operators. One of such determinant representations has 
been obtained in [2] for an equal-time correlator of one-dimensional model of "impenetra- 
ble" bosons, which are described by the quantum non-linear Schrodinger equation with 
infinite coupling. This result has been generalized to the case of correlators with different 
time arguments [3], and also to the case of the XX spin 1/2 Heisenberg chain [1]. The 
determinant representations of the correlation functions allow to deduce the non-linear 
integrable partial differential equations for the correlators pQ, [5]. Various determinant 
representations have been deduced in [U], [7j, [S], 0, [ID1, [Uj (see also Refs. in pQ), 
for instance, for the temperature correlators of quantum non-linear Schrodinger equation 
and for XXX and XX Z spin 1/2 Heisenberg chains. It should be noticed that XX and 
XY Heisenberg models still continue to attract attention [T2], [12], [Ej, and the multiple 
integral representations as well as the determinant representations for the correlators in 
these models are also actively studied [T3], [15]. Multiple integration over a set of the 
Grassmann coherent states is used in ^T], 

In its turn, functional integration (or path integration) technique can be used to cal- 
culate the correlation functions in various quantum models [T7j, [TH], [H], EH]; EB> E2- 
The present paper is based on [22], where an approach has been proposed to represent 
the generating functional of correlators of z-components of spins, as well as the partition 
function, in the Heisenberg XX-model by means of functional integrals defined on the 
variables subjected to, so-called, "automorphic" boundary conditions. 

Approach [22] is based on a technical consideration carried out in Ref.j21] (see also 
|25J) which is concerned with the index theory and supersymmetric quantum mechani- 
cal systems. Path integration is used in [21] to evaluate traces of those supersymmetric 
quantum mechanical operators which appear in dealing with various differential geomet- 
ric indices. In this respect, the Ref. [21] follows [TH], [2H], for the usage of holomorphic 
representation of functional integrals to propose an example of path integral defined on 
the trajectories subjected to non-conventional boundary conditions at the ends of the 
segment of imaginary time. 

It has been shown in [22] that the generating functional of static correlators of z- 
components of local spins in the XX Heisenberg magnet can also be represented by means 
of the (Gaussian) functional integrals which are defined on the trajectories depending on 
the imaginary time non- conventionally in the sense of [21]. In other words, the variables 
of the functional integration are multiplied with a certain complex number when the 
imaginary time is shifted by a period, i.e., the variables behave "automorphically" under 
such shifts. More precisely, the path integrals considered in [22] are defined for the set of 
variables a part of which is subjected to the "automorphic" boundary conditions, whereas 
the other part satisfy the standard requirements of the fermion/boson-type. 

The point is that a trace of an operator exponential is evaluated in |23J, and the 
quadratic operator in the exponent is defined only on the first m sites of the model in 
question (m < M, M is the total number of sites). Remind that the XX model considered 
in [22] can equivalently be handled in the representation of free fermions. Eventually, after 
a passage from a multiple integral over the Grassmann coherent states to the continual 
(i.e., functional) one, it turns out to be possible to define the integration domain so that 
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the integration variables are "automorphic" in the imaginary time on the first m sites, 
while they are (anti-)periodic on the other sites. 

The interest to the functional integrals defined on the trajectories with "automorphic" 
dependence on the imaginary time can also be traced back to |2Z], |2E]- An essential 
distinction between the formulations discussed in [22] and in [21], [21! consists in the fact 
that the "automorphic" boundary condition for the segment of the imaginary time turns 
out to be "inhomogeneous" spatially since it is valid only for a part of sites of the chain 
model in question. The dependence of the integration variables on the imaginary time in 
[27] is the same for all sites. 

The given paper continues [23J, and it is concerned with carrying of the approach 
proposed to the XY Heisenberg magnet which is equivalent to quasi- free fermions (i.e., 
its Hamiltonian in the fermionic representation is diagonalized by the Bogoliubov trans- 
formation). It should be noticed that our path integral representations do not imply a 
straightforward implementation of the proposal [21]: a special restoration of invariance of 
the Lagrangian of the model in question under shifts of the imaginary time by a period 
is required. The method of zeta-regularization is used in what follows to handle the de- 
terminants obtained. The generating functional, as well as the partition function of the 
model, are calculated. Certain correlation functions at nonzero temperature are obtained 
explicitly. Thus it is demonstrated that the path integration approach proposed admits 
a considerable simplification and enough transparency for the problem in question. 

From a physical viewpoint, basic ideas of zeta-regularization (^-regularization) have 
been formulated in [29J, [30], [3Tj, [32] • In mathematical literature, usage of ^-regulariza- 
tion is usually traced to 33]. Zeta-regularization turned out to be rather useful in physics 
to calculate, say, the instanton determinants [31], [35], the Casimir energy on manifolds 
[311], [37], as we ^ as the axial and conformal anomalies [38j. One should be referred to 
[SI], [32], [HB], [HH] for exposition of ^-regularization. 

The paper is organized as follows. Section 2 contains outline of the problem and 
basic notations. The representation for the generating functionals of correlators of a^- 
operators (and also for the partition functions) in the form given by a combination of 
the fermionic functional integrals with "automorphic" boundary conditions is obtained in 
Section 3 (a^ implies the Pauli matrix <j z at nth site). Section 4 contains calculation of 
the functional integrals, i.e., obtaining of the answers in the determinant form. The most 
important formulas of £-regularization are given in Section 5. Moreover, the partition 
function of the XY (and, so, of the XX) model is calculated in Section 5 with the use 
of the generalized ^-function in the series form. The generalized ^-function in the form 
of a Mellin transform is defined in Section 6, and it is used to obtain the regularized 
answers in the form of determinants of finite-dimensional matrices which constitute, in 
their turn, the total generating functional. Differentiation of the integrals obtained with 
respect to a parameter is also considered in Section 6, and some concrete correlators are 
calculated. Reductions of the answers for the XY model to those of the XX model are 
verified. Discussion in Section 7 concludes the paper. 

2 Outline of the model and notations 

Let us consider the XY Heisenberg magnet of spin 1/2 |4L)| . [IT] on a periodic chain with 
the total number of sites M (with even M). Let Q(m) to denote an operator of number of 
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quasi-particles on the first m sites of the chain (m < M). We shall calculate an average of 
the operator exponential exp(aQ(m)) over the ground state of the model (our notations, 
though conventional, correspond to j3j, [TU], [TT]). 

G(a, m) = ($ I e aQ(m) I *o> = ^ (e _^ y) J > « e ^ (2-1) 

where Hxy is the Hamiltonian of the XY model, f3 is inverse temperature (/3 = 

and Tr means trace of operator. The vacuum average (2.1) plays the role of a generating 

functional of static correlators of third components of spins of the magnetic chains jlj, 

0, DDI, HH, US]- 

The Hamiltonian of the XY model has the form: 

H = H + 1 H 1 - hS z , (2.2) 

where 

I M 

H = ~oJ2 °n+l + <*n a n+l) > ( 2 - 3 ) 
Z n=l 

i At 1 M 

Hl = - 5 E (^n+1 + <Wl), ^ = - Y, < , (2-4) 
Z n=l Z n=l 

where S* 2 is the total spin, h is an external magnetic field (h > 0, pQ, jlj, ^3]). The 
algebra of the Pauli spin operators on the sites, a", n G {1, . . . , M}, is defined by the 
commutation relations: 

[a%y n ] = 2i8 kn G a ^al a± = (l/2)« ±<), 

where g q/37 is the totally antisymmetric symbol, and the indices a,/3, 7 acquire the "val- 
ues" x, y, 2;. Besides, the periodic boundary conditions are imposed: ct^+m = °ni Real 
parameter 7 characterizes anisotropy: we obtain the XX Heisenberg magnet at 7 = 0. 

Let us use the Jordan-Wigner transformation from the variables a% to the canonical 
fermionic variables Ck, c\: 

fc-i fc-i 

Cfc = exp (Z7T ^] <7~t7+) CT+, C J = a fc eX P (~ i7r 51 °n ) • ( 2 -5) 

n=l n=l 

The variables c^;, 4 are subjected to the anti-commutation relations: 



{c fc , c n } = {c{, 4} = 0, {c fc , 4} = 5, 



kn 



where the brackets { , } imply anti-commutation. As a result of the transformation (2.5), 
the Hamiltonian (2.2)-(2.4) acquires the following form in the fermionic representation 

H = H + P + + H-P-, (2.6) 

1 M 

H± = —z Yl [ c l c k+i + c{ +1 c k + 7(cfc+ic fc + 44+0. 



2*=! 



M 



+hJ24 c k-hM/2. (2.7) 

fe=i 



The Hamiltonians H ± (2.7) look similiar each to other except for the choice of the spatial 
boundary condition for each of them: the superscripts ± are chosen in correspondence 
with the boundary conditions for Ck, c\ in the form: 

cm+i = =FCi, c ] M+l = tc[. (2.8) 

m , 

We have got in new variables: Q(m) = c I c fc is the number operator of quasi-particles 

fc=i 

on first m sites, the total number of quasi-particles is J\f = Q(M), and the projectors P 
in (2.6) are defined conventionally: P ± = (1/2) (1 ± (— 1)^) 02]- Operator J\f commutes 
only with H (2.3) and S z (2.4) but not with E x (2.4). The parity operator (-1)^ 
anti-commutes with c\ and and it commutes with H . 

As a result, we obtain the following representation for G(a,m) (2.1) |llj : 

G(a,m) = {2Z)-\G + F Z+ + G F Z F + - G B Z B ), 

G F Z$ = Tr UWe-^), (2.9) 

G%Z% = Tr ( e ««W(-l)%^) 

(in what follows we omit the subscript XY in the Hamiltonians). For the partition 
function Z we have got: 

Z = (l/2)(Z+ + Zp + Z+-Z B ), 

(2.10) 

Zj = Tr(e-^ ± ), Z± = Tr({-l) N e-P H± ). 

The following observation is discussed in Ref.j21], which is concerned with the index 
theory and supersymmetric quantum mechanics. Let a, a' be some fermionic canonical 
operators. Let us consider W(l)-operator Q$ = exp(i^a^a), which acts on a, a/ as follows: 

QoaQl = e^a, QWQl = e^aL 

When calculating the trace Tr (Q# exp(— f3H)) in the path integral representation (H is 
a Hamiltonian) , it turns out to be rather natural to come to a path integral over the 
variable subjected to the "automorphic" boundary condition: 

«T) = -e w £(T + 0), (2.11) 

where r is imaginary time, r e [0,/3]. Indeed, Eq. (2.11) reminds the definition of an 
automorphic function (automorphic form, jl^j, jSJ): 

9*f = f(9u) = r(g)f(u) , 

where f(u) is an appropriate function (form), g is an elememt of a group of transformations 
acting on the argument u (thus generating an action of g* on /), and r(g) denotes a 
representation of g*. Thus, Eq. (2.11) implies that the integration variable is transformed 
accordingly to a nontrivial representation of U{1) when r is shifted by the period (3. Other 
physical (quantum-statistical, in fact) examples of "automorphic" boundary conditions at 
the ends of the segment [0,j3] 3 r can be found in [2Z], where spin 1/2 and spin 1 chain 
models are studied by the method of functional integration. 
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It is not difficult to note that operator exp(aQ(m)) behaves analogously: 

(2.12) 



e aQ(m) c e -QQ(m) 



-n i 



c n , m < n < M, 



and it is suggestive to use the idea of [21| when calculating G$Zp, G%Z B (2.9). The "au- 
tomorphic" condition arises for all sites in the models considered in [21]. The peculiarity 
due to (2.12) is concerned with m < M, and the "automorphic" condition is expected to 
appear only for a part of sites. 

To conclude the section, let us define the coherent states using the fermionic operators 
c n , cj, which possess the Fock vacuum | 0): 

c n |0) = (0|4 = 0, Vne{l,...,M}, (0 | 0> = 1. 

Namely, we define the states 



x[a 



)) = exp c l x k(a)^ I 0) = exp(c f x(a)) | 0), 

f M ) 

x*(a) |= (0 | exp < x* k (a)ck > = (0 | exp(x*(a)c), 



;2.i3) 



where a is discrete index running from 1 to N, and the shorthand notations are used: 

M , M 

J2 c k Xk = c } x, n dxji = dx, etc. In fact, independent coherent states are defined 
fc=i fe=i 

which are labeled by independent complex- valued Grassmann parameters x%(a), Xk{a). 
The following relations hold for the states (2.13): 

Ck | x{a)) = Xk{a) | x(a)), (x*(a) \ c\ = (x*(a) \ x* k (a), 

(x*(a) | x(a)) = exp(x*(a)x(a)). 



3 The functional integral 

Let us turn to the problem of rewriting G F Zp, G^Z^ (2.9) and Z F , Zg (2.10) in the 
form of functional integrals. For definitness, let us consider G F Zp\ 



G F Z$ = J dzdz*e z * z (z*\e aQ(m) e-P H± \z), (3.1) 



where it is understood that the trace of operator is calculated as the integral over the anti- 
commuting variables |2E], |45j . and the coherent states (z*\, \z) are defined analogously to 
(2-13): 

= (0| exp(z*c), \z) = exp(c^)|0). 

In order to go over to the path integral, let us divide the segment [0, j3] into iV parts of the 
length P/N, and let us represent exp(— as a product of N identical exponentials. 
Inserting iV decompositions of unity between the exponentials, let us transform (3.1) into 



n r n \ 

G F Z F = / dzdz* Y[ dx*(a)dx(a) exp z*z — ^2 x*(a)x(a) J 

J 0=1 V a=l J 

x(z*\e aQ{m) \x(l)}(x*(l)\e-* H± \x(2)} . . . {x*(N)\e-" H± \z), 



(3.2) 
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where | x(a)} and (x*(a) | are defined in (2.13). 

Using the properties of the coherent states we evaluate the following averages: 



[z e 



where 



and 



(m M \ 

e a J2z* k x k (l)+ Yl *l*k(X) , (3-3) 
fc=l fc=m+l / 

(x*(a)|e""^" H± |x(a + 1)) ^ 
exp (x*(a)x(a + 1) — H ± (x*, x \ a) j , (3.4) 



H (x*, x | a) = H (x*, x | a) + r )H 1 {x* ^ x | a) 
fe=l z 

2 M 

Hq(x*,x I a) = -- ^(4(a)a; fe+ i(a + 1) + x^ +1 (a)a; fc (a + 1)), 

Hf(x*,x | a) = -^53(a;fc+i(a)xfc(a + 1) + 4(a):4 +1 (a + 1)). 
z fc=i 



Inserting (3.3), (3.4) into (3.2), we obtain: 

dzdz* Y\ dx*(a)dx(a) exp <^ ^ ^(^fc + e a 2fc(l)) 

a=l ^ k=l 

M 

+ E 4(^ + a; fc (l)) + x*(l)(x(2)-x(l)) + ... + x*(iV)(^-x(iV)) 

fc=m+l 

-?-(H±(x\x I 1) + . . . + irV, 2 | iV))l. (3.5) 



N 

Let us introduce the notations x k (N + 1) and xjjl(O) as follows: x^(A^ + 1) = ^ (Vfc), 
x£(0) = e a z k (for 1 < k < m) or x^(0) = z% (for m < k < M). Further, we impose the 
boundary conditions with respect to the counting parameter a: 

x fc (0) = -e- a x k (N + 1), l<A;<m, 
Xfc (0) = -x fc (A^ + 1), m < k < M, 

and perform the transition — ► oo. As a result, the discrete index a varying from 1 to 
N then becomes a continuous argument r G [0,/3] (imaginary time; see CH])- The 
right-hand side of (3.5) then becomes the integral: 



/ II dx*(r)dx(T) exp £(r)di 



(3.6) 



where £(r) denotes the Lagrangian: 



£(r) = ^.(r)^ -#*(**, z|r), (3.7) 
fc=i 



and the functional variables Xkij) are subjected to the "automorphic" (see (2.11)) condi- 
tions: 

Xkij) = — e~ a Xk( T + P), 1 < k < m, 

(3.8) 

%k\ T ) = ~ x k{ T + p), m < k < M. 

Generally speaking, the fields x* k {r) are independent integration variables. It is convenient 
to subject x* k {r) to a requirement analogous to (3.8) but with e a instead of e~ a . 

The derivation of the representation (3.6)-(3.8) follows j21] strictly, and it does not 
take into account the peculiar character of our problem: the conditions (3.8) characterize 
two independent sets of sites. Therefore, the following circumstance becomes essential, 
which is new in comparison with [21], [2*7j . 

It can be assumed that certain representations of the group of shifts of r by the period 
(3, i.e., r — > r + (3, are defined by conventinal (anti-)periodicity rules Xk{r) = ±£fc(r + (3), 
k G {1, . . . , M}, as well as by the conditions (3.8). The action functional of the model, 

/3 

/ £(t)g?t, in the exponent of (3.6) is well-defined provided the Lagrangian C{r) is invariant 

o 

under the shifts of r. Such invariance takes place for conventinal boundary conditions 
provided C(r) is even in powers of the fields. 

Let us use (3.8) to calculate the variation 5C(r) at m < M: 

8£{t)=C{t + (3)-£{t) = 

1 



{e a - l)(x* m+l {r)x rn (T) + x* m (t)x m +i(t)) 

+ (e- a - l)(a4(r)x m+1 (r) + x* m+1 (t)x m (t)) 



m—l 



+1 E (> 2a - l)%i(r)it(r) + {e~ 2a - l)4(r)x^ +1 (r) 
z k=i 



+^[(e a - l)(s m+ i(r)a; m (r) + x M+1 (t)x m (r)) 



+ (e a - l)[x* m (T)x* m+1 (r) +x* M (r)x* M+1 (r) / 
The origin of SC(t) is clear: the quadratic forms Hq are invariant under the replacement 

Xk — > ±e a Xk, x* k — > ±e~ a x* k , (3.9) 

provided (3.9) is carried out at each site k 6 {1, . . . , M} (in this case (3.9) looks like a 
homogeneous "gauge" transformation), and they are not invariant provided (3.9) is valid 
only for a subset of {1, . . . , M} (a nonhomogeneous transformation). However, the forms 
Hi are not invariant even for a homogeneous transformation (3.9). The condition (3.8) 
implies a nonhomogeneous representation of the shifts r — ► t+(3, and, thus, the invariance 
turns out to be broken for jC(t) (3.7). However, this symmetry can straightforwardly be 
restored as follows: one should replace H ± (x*,x \ r) in the limiting formula (3.6) by 
another form H ± (t) = H ± (x*,x \ r) of the following type (we omit the superscript ^ at 
H): 

~ ~ M 

H(t) = H (t) + 7 #i(t) + h x* k (r)x k (r) - hM/2, (3.10) 

fc=i 
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where 



^o(t) = ~ J2 ' ( x *k( T ) x k+i(r) +4 +1 (r)x fe (r 
z fc=i 

-^(x;(r)x m+1 (r)e^ + x^ +1 (r)x m (r) e -^ 

+^(r)x M+1 (r)e-^ + x^ +1 (r)x M (r)e^) 

-i m— 1 

z fe=i 

j Af-1 

_ 9 Z (^fe+i(T)x fc (r) +Xfc(r)4 +1 (r)) 



■5( 



fc=m+l 

+a;M + i(r)a;Af(r)e-^/' 3 + ^(<r)x^ +1 (T)e^ 



where the prime at implies that k = m is skipped in summation. The Lagrangian C(t) 
takes the form: 

C(t) = J:xI(t)^P--H(x*,x\t). (3.11) 
fc =i dr 

Equation (3.11) differs from an expected expression since it contains the exponential 
factors exp(±ar//3). The exponentials mentioned ensure periodicity of C(t) with respect 
to the imaginary time [22] ■ Notice that in [21] and [2Z] the "automorphicity" conditions 
are the same at each site, and thus the necessity in the "compensating" factors is absent. 

The Lagrangian (3.11) is invariant under r — > r + /3, the integration measure in (3.6) 
is also invariant, and, finally, we obtain: 

/ P \ 

G$Z$ = J J] dx*(T)dx(r)exp J l{r)dr , (3.12) 



rS[0,/3] 



where C{r) is is given by (3.10), (3.11). For G^Z^ the functional integral representation 
has the same form (3.12). But because of the presence of the parity operator (— 1)-^ 
under the trace symbol in (2.9), we get the corresponding boundary condition in another, 
in comparison with (3.8), form: 

xdr) — e~ a x k (r + (3), 1 < k < m, 

(3.13) 

x k (r) = x k (r + /3), m<k<M 

(i.e., two minus are reversed). Substitution of GpZp, G%Z^, represented by the integrals 
of the type of (3.12) which are supplied with the appropriate boundary conditions, into 
(2.9) gives the desired functional integral representation for the generating functional 
(2.1). 

The main statement of the present paper reads that the representations for GpZp 
and G%Z% given by (3.12) together with the conditions (3.8) and (3.13), respectively, 
are well-defined relations which lead to correct expressions for the correlation functions. 
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Actual calculation below is to argue this assertion. We formally consider (3.8) and (3.13) 
as "automorphic" boundary conditions to distinguish them from more conventional rules 
known for fermions at a = i2ir(k + |) for (3.8), or for bosons at a = i2nk for (3.13), 
k G Z. 

4 Calculation of the functional integrals 

In order to proceed with calculation of the integral (3.12), let us pass to the momentum 
representation: 

x fc (r) = ( / 9M)~ 1 / 2 e~ i7r / 4 £ e i{uJT ~ l f T+qk) x p} 1 < k < m, 

p (A 1) 

x k (r) = {[3M)- 1/2 e- in/4 J2 ei{uiT+qk)x P, m < k < M, K ' J 

p 

where summation goes over the formal 2-momenta p = (oo,q): uj implies the Matsubara 
frequencies uj = irT(2n + 1), n E 2Z (we continue to consider G$Zp), and the quasi- 
momenta q take their values in the sets |llj . 



„, f tt(2Z-1) 
X+ = !q = -7T + 



M 

2tt/ 



1 = 1 M 



(4.2) 



1 = 1 M 



Two sets (4.2) are in correspondence with two "spatial" boundary conditions (2.8). As a 
result of the condition (3.8), the Fermi frequencies in the Fourier expansions for the first 
m sites are shifted by a purely imaginary number. However, the summation index in (4.1) 
can obviously taken as to. We use (4.1) in (3.12) and obtain the following representation: 

GpZp = Y\_dx*dx p exp Sp(a), (4.3) 
J P 

am . 1 , v / %uj-e n -T n 



QpiP2 ^ ]| X P2 

~QpiP2 J \ X -P2 

where the notations are defined: 



(4.4) 



£q h COSq, T^q 'y Sin q, QpiP2 ^uiiiV2Qqiq2i 
Qr 



1 sin f(p - q) 



(4.5) 



^ M sin ^ 



2 



All the quasi- momenta take independently their values in X + or X~ . It has to be 
pointed out that expression S F (a) (4.4) should be viewed a bit formally in the case when 
quasi-momenta belong to the set X~. Indeed, when q G X~, the corresponding Fourier 
coefficients x^, q and require a separate consideration at q = or q = n. All the 
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necessary explanations can be found in [12] (see also [11J). We shall assume that when 
q = or q = it, the corresponding value of the argument — q has to be taken equal 
also to zero or to tt, respectively. Extra term am/2 arises in the matrix representation 
Sp(a) (4.4) as a result of an additional assumption about a fermionic character of the 
anti-commutation relations for the coefficients x* x , x P2 . Compulsory character of such 
assumption is due to the fact that the presence of this term plays an important role for 
obtaining the correct answers for the correlation functions. Calculation for G%Z^ (see 
(2.9)) is carried out just analogously to that above, but the presence of (—1)^ under the 
trace symbol results in the boundary conditions (3.13). Eventually, summation in the 
corresponding quadratic form S B (a) goes over the bosonic frequencies u = 2nTn, n £ Zt, 
and G%Z% looks similarly to (4.3), (4.4) provided all the frequencies are replaced. 

Using (3.10)-(3.12) at 7 = 0, we obtain the corresponding expression in the XX-case 
(Eqs.(4.1) are taken without the factors e~ ln ^): 



G$Zp = J Y[dx*dx p eyLp(Sp(a)), 



q±( \ — W- ^ * a ST n * f3Mh 

Dp{Q>) — } ;{ILU — Eg)X p Xp + — QpiP2 X Pl x P2 "I y 5 

P P Pl,P2 

where e q = h — cosg is the band energy of the quasi-particles defined above. When 
m = M, the matrix Q giq2 (4.5) becomes a Kronecker symbol, and therefore Sp(a) (4.6) 
is transformed into a free fermionic action with shifted chemical potential. 

Let us consider GpZp, given by (4.3), (4.4), and let us carry out the Bogoliubov 
transformation in Sp(a) with 2x2 matrix g e [12], |Hj : 

(x*,x_ p ) = (y*,y- p ) go, (4.7) 

where 

ge = exp(-z e SU(2), tan6 q = 6_ q = -6 q . 

Z E q 

Here it is assumed that a a (a takes values 1, 2, 3) without a lower index implies just a 
Pauli matrix but not a spin operator on a site. Eventually, we obtain the action 



ni / \ am 1^ 

s fw = — + ^2^(yp^y-p) 



iuo — E q a' 



Vp 



L P Pl,P2 \ y -P2 / 

where 0< ee 6 Qi (i = 1,2), and E q = (e 2 q + r 2 ,) 1 / 2 . When q e X~, the Bogoliubov 
transformation is carried out by means of (4.7) at q 7^ 0, tt, while at q = 0, it the remark 
made after Eq.(4.5) has to be taken into account. The point is that, accordingly to 
we put at q = 0, n: 

and 

E = h - 1 = e Q , E n = h + 1 = e n . 



11 



In both cases, q G X + and q G X~, the final answer is written as a single relation (4.8), 
where it is assumed that q and — q are understood appropriately. 

Substitution of (4.8) into (4.3) (with the change of the integration measure) leads j2H] 
to the answer: 



G F Z$ = e am / 2 Det V 2 A( 



a 



A piP2 [a) 



iu F - E qi a 3 \5 PlP2 + - ge 1 {Qp lP2 cT 3 j g^ 1 . 



(4.9) 



Calculation for G^Z^ (2.9) is carried out just analogously to that above, and the final 
answer for G^Z^ looks like (4.9) provided all the frequencies are replaced: 



G%Z% 



e am/2 Bet 1/2 



iui B - E q , a 



PU>2 



a 





502 



(4.10) 



The partition function of the XF-model is given by (2.10) where Zp, Z^ are given as: 



Z$ = Tr (e^ H > 



e-^Det 



[iujp — E q )5 pp , 



where (see [TT] ) 



Z| = Tr 



± 

XY 



e'^oDet 



[iu B - E q )5 pp . 



(4.11) 



X; Eg, E q = (e 2 q + 1 2 sm 2 q) 



Let use the integral representation (4.6) to obtain the following formal answers in the 
XX-case 



Gp^Zp 



e Wh/2 Det 
e PMh/2 Det 



a: 



—iiOp + s q )5 P pi — — Qppi 



-iuj B + e q )8, 



a 



PP 



a Qpp' 



(4.12) 



and 



Z ± = e /3Mh/2 Det ^ iup _ £q)§pp/] 



(4.13) 



pp. 



In Eqs.(4.9)-(4.13) we use the notations ojf and ojb to stress the character, e.g., fermionic 
or bosonic, of the corresponding frequencies. The symbol 'Det' denotes determinants of 
infinite-dimensional matrices, while ' det' is reserved for conventional matrices. Deriva- 
tions of the representations (4.12) for GpZp and G^Z^, as well as of (4.13) for Zp and 
Zq, are carried out analogously. It is convenient to denote the matrix operators, which 
appear in (4.9), (4.10), (4.12), as A(a) = A a , while those in (4.11), (4.13) - as A. 



5 Zeta-regularization (Zeta-functions in the series form) 

We shall use £-regularization jHU] in order to assign meaning to the determinants in (4.9)- 
(4.13). We shall begin with the introductory notes. Usually, a generalized ^-function 
is related to an elliptic operator. Precisely, let A be a non-negative elliptic operator of 
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order p > on a compact <i-dimensional smooth manifold. Let its eigen-values X n being 
ennumerated by the multi-index n. The series 



A„^0 



(5-1) 



which is convergent at 9fts > d/p, defines the generalized (^-function of the operator A, 
((s | A). This series defines ((s \ A) as the meromorphic function of the variable s G (F, 
which can be analytically continued to s = 0. The formal relation 



lim 



o ds 



s \ A) = lim 



log A n 



iog( n m 

\A„^0 / 



allows to define a regularized determinant of A as follows: 



log DetA=- lim ^- ((s I ^4) 

s^o ds 



The Riemann ^-function, 



oo 

C00 = £ n " 

n=l 



and the generalized ^-function, 

oo 

C(s,a) = £(n + a) 



3fo > 1, 



a ^0,-1,-2,... 



(5.2) 
(5.3) 

(5.4) 



ra=0 



are meromorphic in s, have a simple pole at s = 1 with residue 1, and possess a con- 
tinuation at s = 46^ (see also |4*7j). These functions can formally be considered as 
particular cases of the series (5.1). Notice that ( (s, |) is the Gurvitz ^-function [46J, and 

CM) = COO- 

Starting with (5.1), one can represent ((s \ A) as a Mellin transform: 



oo 

£(s | A) = — - J t s - l [Ti (e- At ) - dim(ker^)]dt. 



(5.5) 



The integral (5.5) is defined at sufficiently large positive 9fts (precisely, at 9fts > d/p); for 
other 9fts its analytic continuation is required. Equation (5.5) can be related [HE], [HH| to 
the definition of Det .4 by means of the proper time regularization |48j . |4*H] : 



log 



Det„4 
Det An 



Tr 



dt 

T 



(5.6) 



The definitions of log Det A by means of (5.2), (5.5), and by means of (5.6) coincide up 
to an infinite additive constant. 

Now one can pass to calculation of the determinants (4.11). Let us define the following 
series, which can be expresed through ((s,a) (5.4): 



Ay 



2ni 



s \ '2 2vr 



(5.7) 
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2m 



Wg,q£X ± 



(5i 

- E (-^ 

<7GX± 



i — s 



E 

The series (f( s I ^) an d Cb( s I A) should be considered as the generalized ^-functions 
of the diagonal operators A (see (4.11), (4.13)) in the series form (5.1). The analytic 
continuations for £(s, z) are known 



CM = ~-*, C'(0^)=log^^, (5.9) 



and they lead to the following answers: 

d 

s^o ds 



\imJ-(Hs\A)= J2 Ml + e c ^), 

(5.10) 

-\i m (±( S \A) = E log(l-e c ^), 



where c = ±1 due to an arbitrariness when differentiating (— l) s = exp(±i7rs). 

Choosing c = —1, and combining (5.10) with (5.2), one obtains the following relations 
of the Xy-model PH, HQ: 

Z$ = e~^ J] (l + e"^)= J] 2 cosh 

a rp (5-11) 

Z% = e-^ J] (l-e-^«)= II 2sinh^. 

The total partition function should be calculated accordingly to (2.10), the free energy 
is F = — (1//3M) logZ, while arbitrariness in the choice of c does not influence the 
magnetization M z = —dF/dh and the entropy S = —OF/dT. Specifically, one gets in the 
thermodynamic limit |10j : 



— 

J log(2(l + cosh (3E q ))dq. (5.12) 
o 



2n[3 

All the formulas obtained can be reduced at 7 —>■ to those of the XX-model. 

6 Determinants of the operators A(a) 

6.1 The regular izat ion (Zeta-functions in the integral form) 

Thus, in the previous section we have defined (^-functions of the diagonal operators A 
given by (4.11), (4.13) in the series form. Let us now use (5.5) to calculate the regularized 
determinants of the non-diagonal operators A a given by (4.9), (4.10), (4.12). For instance, 
let us now proceed with the calculation of (4.12). 
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Let us begin with the formal integral 



dt, 



(6.1) 



where e and Q imply the matrices in the momentum space, diag{e g } and Q pq (4.5), 
accordingly, while the indices p, q run independently over X + or X~ . Trace 'Tr' in 
(6.1) is considered as a matrix one over the corresponding 2-momenta p = (u,q) which 
label entries of our matrix operators. Convergence of the integral (6.1) at the upper 
bound is respected at sufficiently large h > h c — 1 (h c is the critical magnetic field [3]). 
Regularization of the integral is necessary at the lower bound. 
Let us use the asymptotical relation 



Tr 



3 (w F -£+fQ)t 



t-+0 <PQ, 



where 0o is an infinite constant equal to Tr (S pp r) = / tr S (5 is a unit M x M matix). 
Let us define the function pit): 

\{iuj F -?+%Q)t 



pit) = Tr 



< t < 1 



(6.2) 



and divide the integral (6.1) into two parts. We rewrite (6.1) using (6.2) as follows: 

oo 1 



r(s) 



dt + 



T(s) 



t s ~ 1 p(t)dt + 



e (iuJ F -E+ZQ)t 

1 v ' 

The function p[t) is a formal series in powers of t n , n > 1. Besides, 

1 



sT(s) 



(6.3) 



sT(s) 



l + 7S + o(s), 7 = -^(l), 



(6.4) 



where if)(z) = (d/dz)logT(z). Therefore, (6.3), which is regular at s — > 0, defines an 
analytic continuation of (6.1) at any > 0. It just can be considered as the definition 
of Cf( s I A*) m the right half-plane oi(F 3 s. 

Let us now consider the constant 0o and the coefficients which define p(t). In our 
situation all these coefficients are given by divergent series, but finite (i.e., regularized) 
values can be assigned to them by means of reductions of the series to zeta-functions (5.3), 
(5.4), i.e., to their particular values at special arguments 

First of all, using £(0) = —\ one obtains: 



M 5^1 = Af(2C(0) 



0. 



where we can equivalently replace ^ and 2£(0) + 1 by and 2£ (0, |], respectively 



( "^-regularized measure" of the set is zero). Further, the divergence of the coefficients 
at the powers of t in p(t) is given by the divergent sums J2 nin - It is reasonable to 
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consider such sums as zeros at m — 2k + 1, k e since n m = + |) are odd. If 
m = 2k, ke % + , then 



E ™ r 



1 N TO 



Z=0 



2C 



-2ifc, - 
' 2 



2 (^-l)C(-2^). 



But ((—2k) — at k > 1. It can be concluded that all "^-regularized" coefficients are 
zero for p(t), corresponding to our A a , and, thus, only the first term is relevant in (6.3). 
Let use (6.4) to pass from (6.3) to the relation: 



d 

lim- C ^MJ = -/ 



tr 



<-?+%Q)t 



dt 



t 



dt 



T - / p(*)t -70o, 



(6.5) 



where the Poisson summation formula enables to sum up over ojp. Then, R.H.S. of (6.5) 
takes the form: 



"trE 



1^ ^ ~/3?+aQ\ fe 



fe=l 



-70o 



logdet (5 + e -/^+«0) - / p(t)| - 70o, 



where 5 is a unit M x M matrix. Let us also take into account that Q 2 
e a Q - 6 = (e a - 1)Q. Therefore, 



Q, and so, 



Dct 



G F 



(iujF — e q )5 pp > + — Q pP 



Det [(iujp — s q )S : 



det 



VP 



5 + (e a - 1)Q(5 + e^y 1 



(6.6) 



Additional renormalization of p(t) and <pQ (to zero, in fact) is irrelevant for G F written 
as the ratio of the determinants. However, when m = M, the corresponding operator A a 
becomes diagonal since Q becomes a unit matrix 5. In this case, we can use ^-function in 
the series form (5.1). Transparent adjusting of (5.7)-(5.10) for this case (i.e., for m = M) 
gives the same answer as that given by (6.5) where pit) and 0o are taken zero. Further, 
we proceed similarly to obtain: 



G% = det 

where we introduced the matrix notation 

M F , B (a) = 



5 + M B (a) 



y - m 

5 ± exp((3e) 



(6.7a) 



(6.76) 



The matrices Mf,b(<x) are of the size M x M; we take + or — in the R.H.S. of (6.7b) 
provided the subscript in the L.H.S. is F or B, respectively. 

In deriving (6.6), (6.7), we have restricted ourselves to the case where h > 1. For 
< h < 1, the energy e q is not strictly positive, and a potential problem of convergency of 
the integral (6.1) at the upper bound arises. It can be shown that in this case it is necessary 
to use the integral representation (6.1) to calculate a regularized value of Det l l 2 {A a A* a ), 
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where A* a is a complex conjugated to A a . Consideration of the corresponding complex 
"phase" of the determinant of A a leads to the same final answer as that at h > 1. 

Let us illustrate the last statements by the following calculation. We are gowing to 
calculate Det (icup + A). Firts of all, let us use (5.5) to calculate 

Det 1,2 (uj f + A 2 ) = VBet(iLU F + A) Det {-iuj f + A), 

where A is nonzero real number. Let us start with the representation 



oo 

hi*- 11 " 



r(s) 



We use the estimation: 



Tr 



e -(u,%+X 2 )t 



P 



dt. 



2(7rf)V2' 

where the Poisson summation formula is used to re-express the series 



^exp(-^U) 



P 



2 J nt 



£(-l) fe exp 



P 2 k 2 ' 
At 



, < t < 1, 



and we define the function p(t) (compare with (6.2)) as follows: 

P 



p(t) = Tr 



2(nty/ 2 ' 



Now, the corresponding generalized ^-function of the diagonal matrix operator diag {up + 
A 2 } takes the form: 



r oo 

C(s diag + A 2 }) ee f t^Tr (e-^ +x ^)dt 

1 (s) J 



+ 



J t s - l p{t)dt 



+ 



P 



27TV2 



(6.8) 



Equation (6.8) defines an analytic continuation of ^-function in the right half-plane of 
W 3 s. 

Using (6.8) we proceed with the calculation of the logarithm of the determinant 
DetV^ + A 2 ): 



i^ c (H dia §^ + A2 >) 

P 



E- 



dt 1 



= log(l+e^ |A| ) 



47TV2 



2. dt 



-x 2 t 



log (2 ch 



f 3/2 ^ 

/3|A| 



- 1 



t 2 



t - o p(*)-t + 







t 27TV2 



t 3 /2 27TV2 



(6.9) 
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In right-hand side of (6.9) we can relace |A| by A. Thus we have obtained the regularized 
value for the logarithm of Det 1 ^ 2 (co F + A 2 ). We determine the "phase" of Det (iu F + A) 
as follows: 

Yl arctan -f = ^T Yl lo S T + = ~ ic ^T' c = ±l, (6.10) 

w F A 2 WF X-tuj F 2 



where we take into account the fact that ^1 = 2£(0) + 1 = 0, and we use 



at £ = ±A. Combining (6.9) and (6.10), we obtain, in agreement with the previous 
calculations (5.7)-(5.10), the following answer: 

Det {iu F + A) = 1 + e c/3A . (6.11) 

Therefore, for any sign of the parameter ("energy") A, we can use a freedom in the choice 
of c so that the final result, say, for the eigen-value matrix e = diag {e q } with e q (4.5) 
will be the same at < h < 1 as that for the magnetic field h > 1. Since our matrix 
operators are diagonalizables, the calculation of (6.11) by means of (6.8)-(6.10) justifies 
our manipulations leading to (6.6), (6.7) at any h > 0. For the bosonic frequencies ujb 
the situation is similiar. 

Now, let us turn to the XF-case. We shall define the regularized determinant of the 
matrix operator A(a) (4.9) by means of the generalized zeta-function of this operator 
((s\A(a)) as follows: 

log Det ^(a) = -Inm^C(sKcO), 

where a standard representation of ((s\A(a)) by the Mellin integral (5.5) is meant. Re- 
ferring to the previous calculations concerning the XX-case, we shall write the answers 
as follows: 

am/2 



detKV«) det& { I + ^LVT^l -y V" } , (6.12) 



Det ^ 


iu F — E 


, a 

8 + ~R 

P 


s(<3 


®<r 3 )# 1 




Det V2 


iujp — E 


<g) <7 3 




) detsfi ' 


\j~9[Q 


® (exp(- 


«(T 3 ) 


- 1 1 


[ J + exp(- 


-/3E C 


5 a 3 ) J 



Z| = det2^{/ + exp(^®cr 3 )}. (6.13) 

The determinants of the infinite-dimensional and finite-dimensional matrices are denoted 
in (6.12), (6.13) as 'Det' and 'det', respectively, while the indices M and 2M imply the 
size of the matrices. In (6.12), / denotes a unit 2M x 2M matrix, g denotes a block- 
diagonal 2M x 2Mmatrix with the blocks on its principal diagonal given by the matrices 
g e defined in (4.7); Q and E are M x M matrices with the entries Q qiq2 and E qi 8 qiq2 , 
correspondingly, and I is 2 x 2 unit matrix. Besides, the following relations are used: 

exp(-a^(Q(g>o- 3 )^ -1 ) -1= g\Q ® (exp(-acr 3 ) - I) 
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and det(expaQ) = exp(am) (see notations in (4.5)). The notations for tensor products 
in (6.12), (6.13) are in agreement with the matrix notations in (4.8)-(4.10). 

Let us verify some reductions of the relations (6.12), (6.13). First of all, let us calculate 
Z$ (6.13): 

/~n~ (l + e^)(l+e-^)= n 2cosh^. (6.14) 

/ qeX± q&X± 2 



7 ± 
F 



Further, we shall consider (6.12). When m = M, the matrix Q becomes a unit matrix S, 
and we obtain at a = in (see (5.11), as well as the definitions (2.9) and (2.10)): 



G$ = J] tanh 



P E 1 - 7±/7± 

—— — Zj b IZj f 



(6.15) 



q€X± 



Let us now put m = M, and a is arbitrary. Then, we obtain from (6.12): 

(1 — ch a)I + sh a g28 c 3 



q± = e aM/2 / -Q 



q&X-- 



D aM/2 



n ± (^ 2 f-sh 



/ + exp(-pE q cr 3 



«cos^th^ + sh 2 ?th 2 ^ 



.16) 



where g2# is the matrix of rotation by the angle 26 q , I is a 2 x 2 unit matrix. On the 
other hand, the result of JT] reads: 



G$ = det M (S + K F [ 



a , 



K F (a) = (e a - 1) | diag {l - e i£ Hh ^}. (6.17) 



Direct calculations of G$ by means of (6.17) taken at m — M, and by means of (6.16), 
coincide. Finally, let us take 7 = while a and m both are arbitrary in (6.12). Then we 
obtain the answer (6.6) for the XX follows: 



Gt = det 



1/2 ( exp(«g) o ^ det i /2 



2A1 







5 



2M 



( ? exp(-aQ) - 5 

o + — — 

5 + exp(— pi) 







V 



o 



exp(ctQ) — 5 
S + exp(Pi) J 



det M ( S + 



exp(aQ) — 5 
5 + exp(Pi) 



(6.18) 



where the matrices Q and S are defined above, and i = diag {e q }, e q = h — cosg. 

It is clear from (6.14)-(6.16), (6.18) that the representations (6.12), (6.13) are in 
agreement with the results (6.6), (6.7) [2*3*| . and they reproduce the relations of the XX 
model jl], [TO], [H], correctly. Calculation for G%Z% is carried out just analogously to 
that above. The final answers for G% and Z% are: 



G± am/2 
B — e 



Det ^ 2 


iwb - £ (8) (J 3 + ^ #(Q ®a^g 1 


Det V2 


— E ® cr 3 
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where 



= detl/ 2 (e^) detJ2{/ + M B (a)}, 

Zp = detlS {l-exp((3E <g> a 3 )} , 

fif[g (8) (exp(-aa 3 ) - J)]^- 1 



M F n(a 



(6.19a) 



(6.196) 



(6.19c) 



J±exp(-/3£® <r 3 ) 

In (6.19c) it is also meant that + or — correspond to F or B, accordingly. We see that 
the answers (6.12) and (6.19a) differ only with respect to the signs at exp(±/3E ® a 3 ). 
Analogously, the results in 11 for G F and Gg also differ only with respect to the signs 
in front of exp(j3E q ). 

It is appropriate to notice that the representation 



G + am /2 Det l/2 J + 



a 



g(Q® a 3 ) g- 1 



(3 iojF — E ® a z 



e am/2 exp 



1 f Tr 

2& k 



a 



[3 iuop — E (g) cr 3 



(6.20) 



(which implies regularization of DetA(a) more conventional for quantum field theory) 
leads to the same numerical coefficients at the powers of a, as Eq.(6.12). For the XX- 
case, (6.20) acquires the following form: 



a 



Gt = Det IS + - 



Q 



(3 %ojf — e q 



^ (-!)*-! 

exp < 2^ : lr 

fe=l 



A: 



a 



Q 



-, k ' 



(3 iuof — £ 



(6.21; 



In (6.20) and (6.21), I and 5 imply the corresponding unit operators. 



6.2 Differentiation of the determinants 

Calculation of the correlation functions by means of the generating functional G(oe,m) 
(2.1) is related with differentiations of it over a and m [JJ, [3], [I], [7j, 0, [TU], [H], [To] . 
In fact, in order to calculate correlators of z-components of spins (the operator of third 
component of spin, a^, is defined as <j z at mth site), we differentiate G(a,m) over a at 
a = as follows [JJ, jt]: 

(g n (m)> = Hm — G(a,m). (6.22) 

However, with regard at (6.1), it is suffice to do only a first differentiation. The other 
ones occur as usual differentiations of matrices [oD] . 

The operators in question, A (a), are linear in a: A(a) = A\ + aA 2 . Let us calculate 
the first derivative of DetA(a) using the formal integral (6.1): 



(d/da)DetA(a) 
DetA(a) 



4(fTr(e^ 
act Wo t , 



(6.23) 
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(the regularization at t \ is irrelevant for the differentiation over the parameter). In the 
spirit of the Ray-Singer-Schwarz lemma |5T] . we shall use in (6.23) the following relation: 

|fTr(e A ^))=4TrfB(a)e^), B(a) = A 2 A~\a). 
act v ' at v ' 

Then, the integral over t can be calculated, and we obtains in the XX-case: 

(1°S Det A a ) — Tr B a 

= Tr ^( iu , F -£+~Q)- 1 \ =tr(Q(l + e^)- 1 ), (6.24) 

where the 'fermionic' A a (4.12) is meant, and the Cauchy formula for matrices (SHI is used 
to sum up over the frequencies. In the Xl^-case, for the matrix operator, say, A a (4.9), 
the answer reads: 

A( logD et^) = -Itr L 9±Q®° 3 )~9-\ \ (625) 

Knowing (6.24), (6.25), one can carry out all the differentiations required. Thus, we 
obtain, 



a) in the XX-case: 



(d/da)Det A a 
Det A a 

(d 2 / da 2 )Det A c 



(6.26) 



lim = tr (Q(5 + e*)-*U 

DetA* ^ v 7 

+tr 2 (Q(o + e^)- 1 ) - tr (q{6 + e^)- l Q(5 + e^y 1 ) , 



etc.; 

b) in the XF-case: 



(d/da)Det l / 2 A a 1 / 5 (Q ® ^ 3 ) _1 



lim — — — — = — tr 



«-o DetV2A a 2 I J + exp(-/?E(g)(T 3 ) 



(6.27) 



,. {d 2 /da 2 )Bet l ' 2 A a 1 / 9 [Q ® I) 9' 1 
hm — — 77— = - tr 1 



Det 1 ^ 2 W + exp(-/?£® a 3 ) 



9 f<3 ® g' 1 



+ T tr ' . 

I \ I + exp(- f3E®a^ 



tr 



(<2 



(q 



^ 3 ) r 1 



7 + exp (- (3E <g> a 3 



I + exp ( — (3E ® (X s 



etc. 
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Equations (6.26), (6.27) are obtained with the help of (6.24), (6.25), i.e., from the 
integral representations of the type (6.1). As an additional check of consistency between 
(6.24), (6.25), on the one hand, and the results given by finite-dimensional determinants 
deduced in Subsection 6.1, on the other, one can verify that appropriate differentiations 
of det(S + Mjr(a)) (6.6) (with Mj?(a) (6.7b)) over a also result in (6.26). Similarly, 
differentiations of det 1//2 (J + M F (a)) (6.12), where M F (a) is given by (6.19c), also lead 
us to (6.27). 

As an other example, let us calculate, in the thermodynamic limit, the derivatives 
(Gp)a and (Gp)" at a = using, for a comparison, (6.12) and (6.17). To proceed, we 
adopt the following definitions: 



C„ 



cos 6 q tYi 



S q = sin q th^ 1 . 



It is straightforward to establish the following relation: 

d 



da 



tr 



a=0 



K(o)) 



m + tr (Mp(0 



m r dq 
' ' G «2V 



where C q is given by (6.28). Further, 
d 2 



da 2 G±F 



a=0 



tr 2 (i^(0)) +tr (k£(0)) - tr (k f (0)K f (0)) 



1 r 
4 



m 



+ tr (M;(0))] 2 + -tr (M;(0)) - -tr (m f (0) M f (0 



We take into account that tr (M/(0)) = m, tr (i^'( )) = tr (-f^(O)), and 

7T 

-Itr(M;(0)M;(0)) = -^ /cos^th^g 



(6.28) 



(6.29) 



(6.30a) 
(6.306) 



^ 7V 7T 

+ 4 / / QpiQ<]p( S p Sc i ~ C p C i) dpdq, 



where C q and S q are given by (6.28). Then, direct use of explicit form of K F (0) to express 
tr (k f (0) K F (0)*) in (6.30a) demonstrates us that (6.30a) and (6.30b) coincide. 

To conclude the section, let us use the obtained formulas (6.29), (6.30) to calculate the 
correlators (cr^) and (<7^ +1 Ci) in the thermodynamic limit. In this case, the contributions 
G%Z~b and G^Z^ in (2.9) cancel and only derivatives of G F are important. We obtain 
for (Q n (m)): 



7T 7T TV 

(Q(™)) = £ / Cqdq, (Q\m)) = £ / Cqdq + / Cqdq 



2 
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+ 



7T 7T 



(6.31; 



where C g , are given by (6.28). The usage of (6.20) instead of (6.12) also leads to (6.31). 
The relations obtained lead to the correct answers for the correlation functions jl], |10j . 
Ej: 



«> = (l/27r) J cos9 q th^dq 



where the definitions 



«) = 1 - 2P 1 (Q(m)) 
and V 1 f(m) = f(m) — f(m — 1) are used. We obtain further: 

7T 7T 

«+i<t1) ~ «) 2 = (l/4vr 2 ) J J cosm(p - q){S p S q - C p C q ) dpdq. 
where the definitions 

(a z m+1 a z ) = 2V 2 (Q 2 (m)) + 2a z -l 

and D2f(m) = f(m + 1) — 2/(m) + f(m — 1) are used. 
In the XX-case we obtain, say, from (6.26): 



(Q(m)> 



m 
2^ 



7T 

J (l + e^^dq, 



(6.32) 



(6.33) 



(6.34) 



TV i 71 \ " 

(Q 2 M) = ^ /(l + e^Y l dq + £/(! + e*-)" 1 *? - 



Q 2 pq {l + e^r'il + e^^dpdq. 
One obtains from (6.34) in the thermodynamic limit: 

z — i z \ i 1 J d( i 



(6.35) 



(6.36) 



The result (6.36) agrees with the magnetization M z = —dF/dh, which is calculated by 
means of F (5.12). We obtain from (6.35) (m > 0): 



z\2 



= if*) 



7T Z 



1 + eP £ 



dq 



(6.37) 



Equations (6.32), (6.33) reproduce (6.36), (6.37) in the limit 7 — > 0, i.e., at 9 q — > 
in C q , S q . The answers obtained (6.32), (6.33) agree with the results of jl], [TO] , 
jllj . [To] . Therefore they witness in favour of self-consistency of the functional- integral 
representations, given by (3.12) (to be used in (2.9)) with the corresponding boundary 
conditions, which appear when calculating, in the way inspired by [24J, the generating 
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functional G(a,m) (2.1) for the XY and XX Heisenberg cyclic chains. The boundary 
conditions in the imaginary time which follow (3.12) are given by (3.8) (for G%Zf) and by 
(3.13) (for G%Z%). Besides, our regularizations resulting in the general expressions (6.6), 
(6.7), (6.12), (6.19) look satisfactory since final results agree with the results obtained by 
other methods. 

7 Discussion 

The generating functionals of static correlators of third components of local spins for 
the XY and XX Heisenberg spin 1/2 chains are calculated in the present paper. The 
calculation is carried out by means of the functional integration over trajectories with "au- 
tomorphic" dependence on the imaginary time. The results are obtained in the form of de- 
terminants of the matrix operators which are regularized by means of zeta-regularization. 
It is demonstrated that, from a practical standpoint (i.e., if only differentiations over 
the parameter a are needed), the formula for the first derivative over a of the generat- 
ing functional can be obtained without regularization of the Mellin integral at the lower 
bound. Various special limits are considered for the formulas obtained which demonstrate 
reductions of the XV-case to the known relations of the XX-case. 

The given paper continues [23J where a method has been proposed which allows to 
calculate vacuum average of an exponential of quadratic operator as a functional integral 
over "automorphic" trajectories. The present paper is close to [21] where the functional 
integration with "automorphic" boundary conditions is used for a calculation of certain 
differential geometric indices. The present paper is also close to |2Z|, where the partition 
functions of spin 1/2 and spin 1 chain models have been also obtained in the form of path 
integrals over variables subjected to "automorphic" boundary conditions. The distinction 
of the present paper from |2*7j consists in the fact the "automorphic" boundary con- 
dition in the imaginary time appears only for a part of sites. The method proposed in 
|23~] is considered above for the system which is equivalent to quasi- free fermions (i.e., the 
corresponding Hamiltonian is diagonalized by the Bogoliubov transformation). Approach 
presented provides further development of the technical finds discussed in j21], j^Z], [2~%] . 
and it seemingly merits attention since can be used further, for instance, for the XX 
Heisenberg model with translationally inhomogeneous boundary conditions. In general, 
the functional integral considered merits attention since it can be useful, as a techni- 
cal method, for other models where it is also necessary to calculate vacuum averages of 
exponentials of quadratic operators of the type of exp(aQ(m)). 
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NOTE ADDED IN PROOF: 

Equation (6.18) demonstrates, in fact, that the square root in G$ (6.12) can be calculated 
at 7 = thus leading to an expression in the form of the determinant of the matrix of the 
size M x M. The same is true for 7 ^ also: calculation of the square root in (6.12) at 
7 7^ leads to G F just in the form (6.17). The corresponding details should be presented 
elsewhere. 
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